A Bose-Einstein condensate bouncing off a rough mirror 
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We present experimental results and theoretical analysis of the diffuse reflection of a Bose-Einstein 
condensate from a rough mirror. The mirror is produced by a blue-detuned evanescent wave sup- 
ported by a dielectric substrate. The results are carefully analysed via a comparison with a numerical 
simulation. The scattering is clearly anisotropic, more pronounced in the direction of the evanescent 
wave surface propagation, as predicted theoretically. 

PACS numbers: 



INTRODUCTION 



EXPERIMENT 



Since the first experimental realisation of a Bose- 
Einstein condensation on a micro-chip pj, the study of 
the interactions between ultra cold atoms and surfaces 
is of crucial interest. For instance, the quality of the 
wires used in micro fabricated chips is directly linked to 
the fragmentation effects observed in trapped BEC near 
a metallic wire Q. Also, it has been shown than the 
thermal fluctuations of the current in a metallic surface 
could induce spin flip losses in an atomic cloud when the 
distance to the surface is smaller than 10 /im typically 

S0. 

The use of dielectric surfaces and evanescent waves for 
producing strong confinement has not been so widely ex- 
plored yet. It has the advandage of a strong suppression 
of the spin flip loss mechanism 0. With such a system, 
one can realize mirrors H , diffraction gratings [g , 2D 
traps |7j or waveguides |8j. Experiments involving ultra 
cold atoms from a Bose-Einstein condensate at the vicin- 
ity of a dielectric surface started only recently, leading 
for instance to the realisation of a two dimensional BEC 
or to sensitive measurements of adsorbate electric po- 
larization jTo| . 

In this paper, we present recent experimental results 
of Bose condensed atoms interacting with the light field 
of an evanescent wave above a dielectric slab. An opti- 
cal waveguide on a dielectric substrate was designed for 
studying two dimensional geometries When a blue 
detuned light is injected into this waveguide, it results in 
an atomic mirror for the condensed atoms released from a 
Ioffe Pritchard magnetic trap. We observe a strong scat- 
tering of the atomic cloud due to the surface roughness. 
Section B is devoted to the description of the experiment 
and the analyse of the results. Section D presents the com- 
parison with a theoretical model developped by Henkel 
et al. [13]. 



The evanescent wave mirror 




Ti0 2 n = 2,3 360 nm 
Si0 2 n = 1 ,46 400 nm 
LaSFN15 n = 1,86 substrate 



FIG. 1: Dielectric prism supporting the evanescent wave. The 
surface is coated by two layers of successively low and high re- 
fraction index to realize a wave guide and enhance the evanes- 
cent field. For each incident polarization, TE or TM, coupling 
is resonant for a given incident angle. The experiments were 
performed with TE polarization. One denotes x as the prop- 
agation axis of the evanescent wave along the surface, y as 
the other horizontal axis and z as the vertical one. 

The atomic mirror is based on the short range repul- 
sion of a blue detuned evanescent wave [f| . The evanes- 
cent wave is produced from a dielectric Ti02 waveguide 
of index n g = 2.3 and thickness 360 nm, separated from 
a glass substrate of index 1.86 by a SiC-2 gap of index 
1.5 and thickness 400 nm (figure^. This surface coating 
was designed to enhance the evanescent field above the 
waveguide 13] . The incident angle is fixed by the reso- 
nance condition for a TE polarization inside the waveg- 
uide. At the Ti02 - vacuum interface, the resulting angle 
is 0i = 46.1° which leads to a decay length of the light 
field K _1 = 93.8 nm where / = I e~ 2tiZ . In the hori- 
zontal plane, the mode propagates along x with a wave 
vector K ev e x with K ev — 1.66fcj,, — 2ir/\ being the 
free space wave vector and Al = 780 nm the laser wave- 
length. 

The light coupled into the waveguide is produced by 
a Gaussian shaped laser diode of power 40 mW detuned 

5P 3/2 ,F' = 3 



1.5 GHz to the blue of the 5S 1/2 , F = 2 
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D2 line of Rb. The Gaussian beam is elliptical and is 
focused at the prism surface to a 1/yfe waist diameter 
of 100 /jm in the plane of incidence and 85 ^m along 
y. Here, one denotes x as the propagation axis of the 
evanescent wave along the surface, y as the other hori- 
zontal axis and z as the vertical one (see figure QJ. After 
projection onto the surface, the x waist is w x — 220 /im 
whereas the y waist remains unchanged. The light inten- 
sity at the center of the spot at the surface is 210 W/cm 2 . 
This value is deduced from a measurement of the light 
intensity threshold for the atomic reflection, taking into 
account the van der Waals attraction to the surface. Note 
that the number of photons scattered by each atom dur- 
ing reflection is about n sp ~ 0.13 only. 

Bouncing off the mirror 

Experimental details on the Bose-Einstein condensa- 
tion set-up are given in Ref. 87 Rb atoms are pre- 
pared and condensed in the F = 2,mf = 2 state in a 
magnetic Ioffe-Pritchard type trap, elongated along x, 
with oscillation frequencies respectively u x /2tt = 21 Hz 
along x and lo±/2tt = 220 Hz in the radial directions. 
The center of the magnetic trap sits at ho = 3.6 mm 
above the evanescent mirror. This distance can be ad- 
justed with of another pair of magnetic coils Q. About 
N = 3 x 10 5 atoms cooled to below the condensation 
threshold are released from the magnetic trap. They ex- 
perience a free fall of 27 ms before reaching the evanes- 
cent wave and being reflected. The laser beam creating 
the mirror is switched on for only 2.2 ms to avoid photon 
scattering during free fall and after reflection. This time 
is nevertheless sufficient to ensure a full reflection of the 
whole atomic cloud. An absorption picture is taken at 
the end of the experimental sequence, either before or 
after reflection. 

From the pictures taken before reflection, we infer ex- 
perimental parameters which will be useful for further 
analysis. One calibrates the pixel size by the center of 
mass motion under gravitation. We have also noticed 
an initial horizontal velocity v x = —30.67 mm/s, due 
to a transient inhomogeneous magnetic field during the 
switching off procedure. From the balistic expansion, 
we get the fraction of condensed atoms (N /N = 0.4), 
the temperature of the thermal cloud T = 285 nK, the 
Thomas- Fermi radius along x R x = 90 /im and the veloc- 
ity radius along z V± = 5.96 mm/s. The velocity radius 
along x, V x , is not measurable, being very small. How- 
ever, following Ref. 0], it can be deduced from V± and 
the ratio of the oscillation frequencies in the magnetic 
trap: V x = § %*-V± = 0.89 mm/s. 

The picture taken after reflection show a strong atomic 
scattering (figure |2fr). The atoms are spread onto a 
sphere with a radius increasing approximately linearly 
with time. From the Gaussian 1/y/e radius of the 



atomic cloud along x, we infer a velocity width of a Vx = 
39.4 mm/s (Gaussian 1/y/e radius), corresponding to 
6.6 v xec where u rec = 5.89 mm.s -1 is the recoil veloc- 
ity for rubidium. Note that the initial velocity width 
contributes to less than 1 per cent to this value. This 
strong atom scattering has no relation with the weak 
spontaneous photon scattering, responsible for a slight 
blur of the picture (on the order of 1 pixel rms of the 
final 512 x 512 image). 




FIG. 2: (a) Absorption imaging pictures of a BEC with 3 x 
10 5 atoms bouncing off a diffuse mirror, for different delays 
after reflection: 2 ms, 7 ms, 12 ms, 17 ms and 22 ms (the 
pictures are superimposed). The prism surface, slightly tilted 
from the imaging axis, is visible at the bottom of the picture. 
(b) Simulation of a bouncing BEC in the same conditions, 
with a velocity spread along y chosen to be a Vy = a Vns /2 = 
19.5 mm.s -1 . The position of the mirror surface is marked by 
a grey line. The dimensions of both pictures are 5.7 mm x 
4.4 mm. 



Analysis of y spreading 

It is more difficult to have access to the velocity spread- 
ing a v along y, as it corresponds to the imaging axis. 
Nevertheless, it can be deduced from a careful analysis 
of the data, in conjunction with a numerical simulation 
reproducing the experiment for different possible values 
of the y velocity spreading. It is clear for instance that, 
in the case of a totally isotropic scattering, the atomic 
cloud would shade on the xz plane towards the surface. 
In contrast, a one dimensional scattering along x would 
present a circle line. 



Numerical simulation 

The experimental picture is compared with a picture 
calculated for a diffuse reflection, knowing all initial pa- 
rameters of the expansion and the amplitude a Vx of the 
scattering along x. The corresponding y amplitude a v 
is given as an input parameter and the output image 
should correspond at best to the experimental image. 
N = 3 x 10 5 classical atomic trajectories are calculated, 
with random initial positions and velocities which match 
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the initial Thomas-Fermi position and velocity distribu- 
tion for 40 % of atoms in the condensate, the other dis- 
tributions corresponding to those of a thermal cloud at 
285 nK in the harmonic trap with known oscillation fre- 
quencies. Only the initial position of the condensed part 
width along y and z is neglected, because its contribution 
to the total width is very small already after a few mil- 
liseconds of ballistic expansion. In the calculation, the 
mirror is flat and the z component of the atomic velocity 
is instantaneously reversed, thus neglecting the exponen- 
tial character of the real mirror: it doesn't play an impor- 
tant role in the dynamics at a such height, the bouncing 
time being in the order of the microsecond. The mir- 
ror roughness is taken into account by adding a random 
horizontal velocity to the incoming one, distributed on a 
Gaussian with a fixed width a Va . along x and a tentative 
width a v along y. The vertical velocity of each atom is 
modified correspondingly to ensure energy conservation. 
The final calculated image is produced by integrating the 
atomic density along the imaging direction (y axis). 

Deducing the y spreading 

The direct comparison of the experimental and calcu- 
lated picture confirms that the simulation nicely repro- 
duces the overall behaviour of the atoms (figure ■ To 
be more quantitative, we isolate a small region of size 
0.8 mm x 1.5 mm along x and z respectively, centered 
on the peak density, and compare the profiles deduced 
after integration along x for both pictures (the integra- 
tion cannot be performed over the whole sample due to 
the curvature of the signal). Figure |31 gives an example 
of such profiles, after a 59 ms total time of flight. 

As can be seen on the figure, the experimental results 
are not consistent with an isotropic scattering, neither 
with a fully ID scattering. A reasonable value for tj v 
would be cr Vx /2 = 3.3 i> rcc , that is an anisotropy x = 
o« / a v x = 2. We cannot exclude values of x as l° w as 
1.5 or as high as 2.5; in any case, the scattering is larger 
in the direction of propagation of the real part if ev e x of 
the optical wave vector. 

THEORY 
Presentation of theory 

The problem of diffuse reflection of atoms off an 
evanescent wave mirror has been addressed theoretically 
in Ref 01 an d it linked to the scattering of light by the 
defects of a rough dielectric surface. Atom scattering is 
due to Raman stimulated transitions involving a photon 
of the evanescent mode and a scattered photon. The 
expected momentum distribution after reflection can be 
calculated from the knowledge of the roughness power 
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FIG. 3: Comparison of experimental results and simulation, 
after 59 ms total time of flight: atomic density profiles inte- 
grated along y and averaged over 1.5 mm along x. Closed 
circles: normalized experimental data. Lines: normalized nu- 
merical calculation, for different tentative values of the scat- 
tering amplitude along y a Vy . The parameters for the simu- 
lation coincide with the measured experimental parameters: 
3.59 mm initial height above the mirror, 40 % of condensed 
atoms, Thomas- Fermi radii V y = V z = V± = 5.96 mm.s , 
V x = 0.89 mm.s" 1 , R x = 90 fim, temperature T = 285 nK, 
initial velocity v x — —30.67 mm.s -1 , v z = 0.3 mm.s -1 and 
scattering along x tr Vx = 39 mm.s -1 . Thin line: a v = (ID 
scattering); dashed line: a Vy = <j v ^/2 (medium anisotropic 
scattering); bold line: cr v = a Vlc (isotropic scattering). 



spectrum -Ps(Q) of the surface and an atomic response 
function B a t(Q) [Eqs.(5.6, 5.7) of Ref.[I3]. Q is an in 
plane two dimensional wave vector. For normal incidence 
onto the surface, the variance of the velocity distribution 
after the bounce is then given by 

^T = TI /d 2 QQ^s(Q)|Sat(Q)| 2 , i = x,y. (1) 

v rec K L J 

J\f is a normalization factor, related to the scattering 
probability w through 

M = w/ Jd 2 QP s (Q)\B at (Q)\ 2 (2) 

Equation JIJ allows one to compute the anisotropy of 
the velocity width after reflection, the dependence on 
the angle of incidence being actually negligible for our 
parameters (see Ref.[l^|h 

The roughness power spectrum may be inferred from 
a measurement of the surface. For most surfaces, this 
spectrum follows a power law -Ps(Q) oc 1/Q a , with 
a being typically between 2 and 5. The atomic re- 
sponse function -B a t(Q) has a complicated expression. 
Nevertheless, for atomic incidence normal to the sur- 
face, it is very peaked around two circles centered at 
±if ev e x of radius Ul- We will thus model this func- 
tion by a delta function around these circles, |B at (Q)| 2 cx 
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S(\Q - K cv e x \ - k L ) + S(\Q + K cv e x \ - k L ). Note that 
this corresponds to the absorption of a photon from the 
guided mode and the successive stimulated emission into 
a free mode at grazing incidence - or vice-versa. These 
two processes contributes to the same amount to cr^., as 
Ps(Q) on ly depends on the modulus of the transferred 
wave vector. The integration may be performed around 
the value Q = if ev e x , leading to: 



4^ - [d 2 Q (Q x - K CV ) 2 P S (Q - K cv e x )S(Q 



k L ) 
(3) 



^ ~ ^ fd 2 Q Q 2 y P s (Q - K cv e x )S(Q - k L ) (4) 

w rcc fc L J 

After integration and simplification, with a power law for 
the roughness spectrum, we get: 



2 _ ^v*_ _ JQ 

X o 

erf, 



dp 



dip 



(cos ip — T]) 2 



(1 - 2-qcosip + rj 2 ) a / 2 

'■ 2 

sm ip 
(1 - 2r] cos ip + r\ 2 ) a l 2 



(5) 



where rj = K cv /kL = n g sin 9i. This expression depends 
only slightly on the choice of a: \ lies between 2.05 and 
2.31 for values of a between 2 and 5 which are relevant 
for a typical polished dielectric surface. It is minimum 
for a — 3 and has the same value for a and 6 — a. In 
the case a = 4, the expression for \ 2 is very simple as 
we have x 2 = 1 — ^V 2 > which for our parameters leads to 
X = 2.12, in good agreement with the experimental data. 



Estimation of momentum spreading 

The evaluation of the expected momentum spreading 
is more delicate, as a full knowledge of B at (Q) and Ps(Q) 
would be required to estimate w. However, if one knows 
the rms surface roughness, it is possible to give an upper 
bound of the effect using the upper limit for the total dif- 
fuse reflection probability w max (Eq. (5.21) of Ref . (l2^1: 



w < w„ 



(4w s e KZ 7A dB r 



(6) 



Here AdB = h/mv z is the mean de Broglie wave length 
of the incoming cloud, z$ is the height of the classical 
reflection turning point, and a s is the surface roughness 
defined as the square root of the variance of the surface 
height. The velocity spreading may then be estimated as 
being, for example in the x direction: 



(T,„ = v r ^w 



(cos ip — rf) 2 
o ^(1 - 2-qcosip + if) a / 2 



rcc "'max r ir 



dip 



A surface roughness of a s = 3.3 nm was deduced from 
an AFM measurement of a portion of the prism. For 
a = 4, and a falling height of 3.6 mm, we find an expected 
maximum velocity spreading along x of 8.00 v rC c, a bit 
higher than the observed 6.6 u rec - 



CONCLUSION 

In conclusion, we have performed measurements on 
the diffuse reflection of a Bose-Einstein condensate on a 
rough surface. The experiment shows a strong scattering 
of the atomic cloud, due to elastic Raman photon scat- 
tering involving both the mirror light and scattered light. 
A first clear evidence of the anisotropy of this diffuse re- 
flection is given 0|. The scattering is twice along the 
direction of propagation of the evanescent mode, and this 
value is in good agreement with the prediction of Henkel 
et al. 12], in a simplified model. A tentative absolute 
value for the scattering amplitude was also given, in a 
reasonable agreement with the experiment. However, a 
more precise analysis would be necessary to confirm the 
anisotropy value and the absolute value 0] . 

The experiments were performed with a Bose-Einstein 
condensate. However, the coherence of the source is not 
visible on the results - no speckle is visible for instance. 
An estimation of the size of the speckle spots, on the 
order of AdB^o/-R where ho = 3.6 mm is the distance 
between cloud and surface and R ~ 100/J,m its size at 
the reflection point, would give typically 1 /im, which is 
not resolved by the imaging system. Nevertheless, the 
reflection should be coherent, the spontaneous photon 
scattering being negligible. To prevent diffuse reflection, 
a great care should be taken in the surface quality. We 
point out however that the atom scattering did not mask 
the diffraction pattern in a temporal diffraction experi- 
ment we performed with the very same prism [lsj . Inter- 
ferometry with evanescent wave devices remains a very 
promising alternative. 
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